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1 Introduction 

The aim of this paper is to apply the theory of spectral curves to analyze the phase structure 
and the critical processes arising in N = 1 SUSY U(N) gauge theories with adjoint matter $ 
obtained by deforming the TV = 2 theories by a tree-level superpotential Tr W(<&). Spectral 
curves with s cuts are associated to the classical vacua that break the gauge group U (N) as 
a direct product of s factors U (N\) X • • • X U (N s ) where N = N\ + • • • + N s . These spectral 
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curves arise as a consequence of the large N dualities between supersymmetric Yang-Mills 
theories and string models on local Calabi-Yau manifolds X of the form [1-4] 

W'(z) 2 + f(z) + u 2 + v 2 + w 2 = 0, (1.1) 

where W(.z) and f(z) are polynomials 

W(z) = ^—+t n z n + ---+t x z, (1.2) 
n + 1 

f(z) = b n ^ l z n - 1 + --- + b . (1.3) 

The corresponding tree-level prepotential T can be characterized as a function of the partial 
't Hooft parameters S% by the special geometry relations 

Si= i n, f£ = I n, (i.4) 

where £1 is the holomorphic (3, 0) form in X and Aj and Bj form a symplectic basis of three- 
cycles in X. In turn, these three-cycles can be understood as fibrations of two-spheres over 
paths in the Riemann surface defined by the spectral curve 

y 2 = W'(z) 2 + f(z). (1.5) 

Integration of over the fibers reduces the integrals (1.4) to integrals of y(z)dz over the 
projections of the fibers onto the spectral curve. 

In this paper we consider spectral curves £ determined from the following data: 

1. A set of s pairs of branch points af joined by s finite disjoint cuts 7$. 

2. A set of s nonzero complex numbers Si such that the branch of y(z) with asymptotic 
behavior 

y(z) = W'(z) + 0{z- 1 ), z^oo, (1.6) 
satisfies the period conditions 

y(z)dz = -AmSj, (1.7) 

where Aj is a counterclockwise contour encircling the cut jj. 

When we need to make explicit the fact that the spectral curves depend on both the cuts 
and the partial 't Hooft parameters we denote these spectral curves by £(7, S), where 
7 = 71 U 72 U . . . U 7 S and S = (Si, S 8 ). 

We emphasize that our analysis of spectral curves does not rely on random models of 
matrices with eigenvalues constrained to lie on some path T in the complex plane (holomor- 
phic matrix models). Our motivation is that a proper definition of these models and their 
planar limit involves several deep subtleties [5]. For example, the saddle point solutions 
which provide the planar limit of the free energy exist only if the path T is such that the 
corresponding eigenvalue density is real and positive. This condition is trivially satisfied by 
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hermitian models, where T is the real line and the coefficients of PF(^) are real numbers, 
but it represents a quite non trivial requirement for general holomorphic matrix models. 
As a consequence, the characterization of the cuts of the spectral curve in terms of the 
support of the eigenvalue density of holomorphic matrix models is a difficult problem. But 
the precise form of the cuts is obviously required to analyze critical processes such as cut 
splitting [4, 6-10] and to study global features of the phase structure of the set of spectral 
curves. 

In this paper we avoid any ambiguity by consistently using as cuts the projections onto 
the spectral curve of minimal supersymmetric cycles in the Calabi-Yau space X [11-15]. 
These minimal cuts are characterized by the condition that the phase of y(z)dz is constant 
along each cut. Minimal cuts are the natural generalization of the spectral cuts of the 
asymptotic eingenvalue density in hermitian matrix models [5, 16]. 

One of the main reasons to introduce matrix models in the study of gauge/string 
dualities is that the planar limit of the matrix model free energy provides the tree-level 
prepotential T [2]. However we will show that for any spectral curve the same expression 
of the prepotential 

T= [ W{z)p{z)\dz\-\ f \dz\ f \dz'\Log(z-z') 2 p(z)p(z') (1.8) 

J -y ^ J y Jy 

naturally appears from a variational characterization of the complex density p(z) defined 
by 

p(z)\dz\=y(z + )— = -y(z^) — , z£ 7 , (1.9) 
and constrained by the period conditions (1.7), which in terms of p(z) are 

p(z)\dz\ = S j . (1.10) 

Incidentally, the subindices z + and z_ in (1.9) refer to the one-sided limits of the corre- 
sponding function y(z) on 7, and in (1.8) the logarithm Log(z — z') 2 has to be understood 

as 

Log(z - z) 2 = \og{z + - z) + log(z_ - z'), z,z'ej, (1-11) 

for consistently chosen branches of log (2 — z'). The prepotential functional (1.8) is inde- 
pendent of the precise form of the cuts as long as they remain in their respective homology 
classes in the complex plane with all the branch points deleted. 

The fundamental application of the prepotential J 7 as a function of the coefficients t of 
W(z) and the partial 't Hooft parameters S is the determination of the vacuum expectation 
values (vevs) 

u fc = i(Tr$ fc ), Si = (Si), k = l,...,n, i = l,...,s, (1.12) 

in the vacuum states |k) labelled by k = (/cife • • • k s ), (hi = 1, . . . , iVj), which correspond to 
the broken gauge group U(N{) x • • • x U(N S ). Thus [1, 2] if we introduce the super potential 
iy eff (t,S,A 2 ) by 

^ eff = E^|^ + 51ogA 2iV , (1.13) 

i=l 1 
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where 

s 

s = Y, s * ( L14 ) 

and A is the nonperturbative scale in the M = 2 gauge theory, the vevs Sj are the solutions 
of the field equations 

__ = 0, * = l,...,a, (1.15) 

and 

Ufe = ^tT' 2.^ = aiogA^' (L16) 

where 

Wi ow (t ! A 2 ) = W efT (t,5 1 (t,A 2 ),...,«S s (t J A 2 ),A 2 ), (1.17) 

is the low-energy superpotential. The logarithm in (1.13) is only defined modulo 2tti so that 
Si and Wiow are multivalued functions of A 2 . An important problem is the characterization 
of the quantum parameter space _A/f q [6, 17, 18] on which all these functions are single- 
valued. For a fixed tree-level superpotential W(z) of degree n + 1 there is a decomposition 
of the quantum parameter space 

M q = M^U---M^ (1.18) 

(s) 

into sectors Mq corresponding to vacua with a fixed number s of factors of the broken 
gauge group. Each point (S±, . . . ,S S ) £ Mq determines a class of s-cut spectral curves 

with partial 't Hooft parameters equal to Si. The characterization of the subsets of spectral 

(s) 

curves with minimal cuts in the sectors M\ and their possible interpolations by smoothly 
varying the parameters of the theory is also an important issue. 

In this work we use a combination of analytic and numerical methods to study the 
phase structure and the phase transitions in the space of spectral curves with minimal 
cuts for a given polynomial W{z). These phases are labelled by the number s of cuts and 
are described by manifolds with points parametrized by the partial 't Hooft parameters. 
Moreover, using the correspondence between points of the quantum parameter space and 
spectral curves we translate our analysis of the phase structure of spectral curves to the 
quantum parameter space of M = 1 SUSY U(N) gauge theories. 

The layout of this paper is as follows. In section 2 we explain our method to charac- 
terize spectral curves and minimal cuts, briefly review some results concerning the classical 
limit, and make precise the notion of critical spectral curves. Section 3 is devoted to the 
study of the prepotential associated to a spectral curve via a variational characterization of 
the complex density; then we consider the corresponding superpotential and the character- 
ization of quantum vacua in terms of solutions of the field equations. In section 4 we study 
critical spectral curves and, in particular, we analyze the phase transition corresponding 
to the splitting of one cut in the quantum parameter space. Sections 5 and 6 contain our 
numerical and analytic study of the spectral curves and the quantum parameter space for 
the cubic model. Using our analytic condition (derived in section 5.3) satisfied by criti- 
cal spectral curves, we characterize the transition curves between the one-cut and two-cut 
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phases in both the space of spectral curves and the quantum parameter space. The paper 
ends with a brief summary and we defer to two appendixes some technical proofs. 



2 Spectral curves and minimal cuts 



In this section we discuss the characterization of spectral curves (1.5) with minimal cuts. 
We denote by 01, . . . ,a n the critical points of W{z) 



W{z) = (z-a 1 )---(z-a n ). 



(2.1) 



We will assume that the roots of y 2 (z) are either simple or double, and denote the double 
roots by a%, . . . , a r . Hence the function y{z) for an s-cut spectral curve can be written as 



y(z) = h(z)w(z), 



where 



h(z) = Y]_(z - ai), w(z) 



i=i 



\\{z - a m )(z - Clm), 



rn=l 



r + s = n, and where the branch of w(z) is fixed by 

w(z) ~ z s , z — > oo. 
But according to (1.6) the factor h{z) in (2.2) is given by 

'w'(zy 



h(z) 



w(z) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



Therefore the function y(z) is completely determined by its branch points, the simple roots 



+ 



1 , . . . , of . We will often use the variables 



ft 



a- — a- 



(2.6) 



2 2 
2.1 Determination of spectral curves with minimal cuts 

To determine the endpoints af, . . . , af of an s-cut spectral curve we substitute (2.2) into 
the left-hand side of (1.5) and identify the coefficients of 1, z, ■ ■ • , z n+s ~ 1 in both members. 
The remaining coefficients do not give independent relations because 



y(zY-W'{z) 2 



W'(z] 



W'(z) 



w(z) ) e L\ w(z) ) e 
Oiz^ 3 - 1 ), z^oo. 



w(z) 2 -2W'(z)w{z) 



(2.7) 



Thus we find n + s equations which, however, involve the coefficients bo, ... , b n -\ of the 
polynomial /. The s additional independent relations follow from imposing the period 
conditions (1.7). 
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Incidentally, it will be useful later to note that as a consequence of (1.6) and (1.7), the 
coefficient b n -\ of f(z) is given by 

b n -i = -4S. (2.8) 

Note also that implicit in the calculation of the endpoints for s > 1 is the choice of the s 
cuts 7i , . . . , 7 S connecting the unknown pairs of endpoints of , . . . , af . 

In general this method leads to several families of solutions for the set of branch points 
as functions of the 't Hooft parameters S which, in turn, determine several families 
of spectral curves £(7, S). Moreover, the values of the branch points depend only on the 
homology classes of the cuts jj in C \ {af , . . . , af}. 

The next step after determining the branch points is to characterize the set of values 
of S for which there exist minimal cuts jj along which the phases of y(z)dz are constant 
(cf. figure 1), i.e., each cut 7^ must be a Stokes line determined by 

ReGj{z) = 0, z£jj, (2.9) 

where 

Gj(z) = e - iarg ^ / 2 y(z' + )dz'. (2.10) 

JaJ 

Note that although the integration path in (2.10) is left undefined, due to the period 
conditions (1.7) the real parts of the functions Gj{z) are single- valued. 

Unfortunately, except for certain symmetric cases only a few general facts can be stated 
about Stokes complexes [19]: since the endpoints are simple zeros of y 2 (z), three Stokes 
lines stem from each branch point forming equal angles 27r/3, and each of these Stokes lines 
does not make loops and ends either at a different zero of y(z) or at infinity. We anticipate 
that the study of critical configurations requires the calculation of the complete Stokes 
graph, which comprises both the Stokes lines (2.9) stemming from the simple roots and 
the four Stokes lines stemming from each double root on (if they exist). The calculation of 
a Stokes graph in a generic case has to rely on numerical methods. 

2.2 The Gaussian model 

The simplest example of spectral curves with minimal cuts is provided by the Gaussian 
model 

W(z) = ^, (2.11) 

for which only one-cut spectral curves may arise. Substituting y 2 = (z — a~~){z — a + ) and 
f(z) = -45 into (1.5) we get 

a + = -a-=2\/S. (2.12) 
Then, for z on the straight line segment 7 with endpoints 

y{z + ) = i|z 2 -4S| 1/2 e iarg5/2 (2.13) 

and therefore 7 is clearly a minimal cut. 
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Figure 1. Stokes lines ji and contours Ai. 
2.3 The classical limit 

A s-cut family £(7, S) of spectral curves is said to admit a classical limit if as S — > the s 
cuts shrink to s non-degenerate critical points of W(z) 

af-aj, (j = l,...,s), (2.14) 

while the double roots of y 2 {z) tend to the remaining critical points of W(z) 

ai~a s+ i, (Z = l,...,r), (2.15) 

so that the family of spectral curves degenerates into y 2 ~ W'(z) 2 . In this case — aj | -C 
\ a j — Ofc| for all k 7^ j and then (2.14)-(2.15) imply that for z £ 7,- 



y 2 (z) = H(z-a l ) 2 H(z-ar)(z-at) 

1=1 i=l 

W"( aj ) 2 (z-a+)(z-aj). (2.16) 



Then 



W"(a 3 ) 



(a- — a. 



\2 



Hence from the period relations we get 



16 W"{a jt 



(j = l,...,s). (2.18) 
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This means that the solutions of the cut endpoints equations in the classical limit is 

s *~W^7y C7 = l.-,'). (2-19) 

Moreover, by adapting an argument used by Bilal and Metzger [5] in the context of 
holomorphic matrix models, we can prove that in the classical limit these spectral curves 
have minimal cuts which, to first order in S, are the straight line segments with endpoints 
op. According to (2.16) 



where 



w "( a j) a*<_ a 7~ a 



(2.21) 



\W"(ajW \aj-aj\ 

Then setting 2ij)j = arg Sj — ifj 7 the period conditions (1.7) imply that the segments [aj , Oj] 
are minimal cuts. An alternative argument to see this property follows from an observation 
of Felder [16]: as S — > 0, the minimal cuts jj are the level lines 

Re ( e -™*gSiW(z) - e~ iaigS m(aj)) = 0. (2.22) 

Hence, if aj is a non degenerate critical point of W(z), any sufficiently small circle C around 
aj intersects the level lines at four points. For small S the point aj splits into the two branch 
points a^ and, by continuity, they must be connected by a Stokes line (2.9) (i.e., a minimal 
cut) inside the circle C so that there are four Stokes lines leaving C. 



3 Prepotentials, superpotentials and vacua 
3.1 Prepotentials associated to spectral curves 

The purpose of this section is to show how prepotentials can be associated to generic spectral 
curves. To this aim we first introduce a semi-infinite oriented path T containing the cuts as 
shown in figure 2, and such that for each z 1 in V there exists an analytic branch of log (2 — z') 
as a function of z in C minus the semi-infinite arc Y z > of V ending at z' , that verifies 

log(z+ — z ) + log(z_ — z ) = log (2+ — z) + \og{z'_ — z), for all z ^ z' in V. (3.1) 

The property (3.1) is essential for the consistence of the following definition of log(z — z') 2 

log(z - z'f = log(z+ - z') + \og{z- - z'), z, z' e T, (3.2) 

which is assumed in the expression of the prepotential T . It is easy to prove that for arcs 
r with parameterization z = z(t) such that at least one of the functions x(t) and y{r) is 
strictly monotone, the property (3.1) is satisfied by the logarithmic branches defined by 

\og{z-z') = \og\z-z'\+ I — — n z'er, zeC\rv, (3.3) 
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7i 



Figure 2. Example of a semi- infinite path T that contains the cuts 7, and allows the construction 
of consistent determinations of the logarithms in the prepotential. 



where Y zz i is any path in C \ T z i connecting z' + \z — z'\ to z. For example if V is a 
real interval of the form (—00, xq] then (3.3) determines the principal branch of log(z — z') 
and (3.2) gives log(z — z') 2 = 2 log \ z — z'\. 

Using hereafter these logarithmic branches, we consider the function 



9W 



\og{z — z')p(z')\dz'\ 



(3.4) 



where p(z) is the complex density (1.9). Note that both y(z) — W'(z) and —2g'(z) are 
analytic in C \ 7, vanish as z — > 00 and, according to (1.9), have the same jump on 7. 
Therefore 

y(z) = W'(z)-2g'(z), (3.5) 

and using 

y(z + )+y(z^) =0, zG 7 , (3.6) 

we get 

W'(z)-(g'(z + ) + g'(z_)) = 0, z£j. (3.7) 

Equation (3.7) means that H^(z) — (g(z + ) + g(z-)) is constant on each connected piece of 
7 or, equivalently, that there are (not necessarily equal) complex numbers Lj such that 



W(z) - (g(z + ) + 5 0_)) = U, ze 7i . 



(3.8) 



A straightforward calculation using (3.1) and (3.2) shows that (3.8) can be written as 
the variational equation 



S_ 



T + J^LifSi [ dq(z)) 
i=l \ J v* ' _ 



0. 



(3.9) 



where 



W(z)dq{z)-\ [ dq(z) [ dq(z') Log(z-z') 2 , 



is the prepotential functional and where 

dq(z) = p(z)\dz\. 



(3.10) 



(3.11) 



Thus, the variational equation (3.9) characterizes p{z) as a (in general, local) extremal 
density for T constrained by (1.10). It also follows that (cf. Appendix A) 



so that we can express the superpotential in the form 



W e z = Y J N i L l + S\ogk 



2N 



(3.12) 



(3.13) 



Note that if we write (3.10) as 



W(z)- l -{g{z + )+g{z^))\dq(z) 



(3.14) 



and use again (3.8), we obtain the usual alternative expression for the prepotential 

T = \ ! W(z)dq(z) + ( 3 - 15 ) 
3.2 The Gaussian model 

As an illustrative example we consider again the Gaussian model (2.11) with the minimal cut 
7 given by the segment [— 2^/S, 2^fS\. Let us take the path V as the semi-infinite straight 
line containing 7 and ending at 2^fS, and define log(z — z') for z not in Y z i according 
to (3.3). Then we have 



log(z + — z ) + \og(z- — z') = 2 log \z — z'\ + i arg S, z, z G 7. 
Hence if we parameterize 7 by z(t) = 2t\S\ 1 / 2 e iarg5 '/ 2 with —1 < t < 1 we get 

4 |£|2 e i2argS pi 



(3.16) 



T 



[ t 2 VY^dt-i^s 2 

71" J-i 2 
4 | CI 2 „i2ar g 5 T rl rl 

11 J J VT^^i - (t>) 2 log \t - t'\dtdt' 



7T- 



+^iog ( 2|sn 



3 1 



4 2 

where log 5 = log |5| + i arg S. 



logS) S 2 



(3.17) 
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Hence the superpotential is 

W cS = NS(1 -logS) + SlogA 2N . (3.18) 
The field equations reduce to 

log(-^)=0, (3.19) 

and we get the N vacua | k) of the J\f = 1 SUSY U (N) gauge theory which are characterized 
by the vevs : 

S« = C fe A 2 , (3.20) 

and the low-energy superpotentials 

W^(A 2 ) = N( k A 2 , (3.21) 

where 

Ck = e 2 ^ N , k = l,...,N. (3.22) 

Thus the quantum parameter space _A/f q is made of ./V copies of the complex A 2 -plane 
connected at the point A 2 = 0. 

3.3 Vacua in the classical limit 

The leading approximation of Sj in the classical limit can be determined [1]. In fact, it 
follows from the the asymptotic formula (7.4) of Appendix A 

Lj ~ W{a 3 ) -J^Sk log A% + Sj (l + log (^^) ) , (3.23) 

where Aj^ = aj — a,/-. The corresponding classical limit for the superpotential is 
s s W"(n \ N i TT A~2A r iA2Af 

w cS ~ y: N i +E 5 4 ^ ( 3 i 13 — ) + M , (3-24) 

and the field equations in this approximation read 

/W"(aj) N i Yii+j A~ 2Nl k 2N ^ t 

log ( ) = 0. ( 3 - 25 ) 

which give rise to N\ x N2 x • • • x N s vacua |k) = \k\k2 ■ ■ ■ k s ) with approximate associated 
vevs [1] 

5 (k) ^ ^ kj ,N 3W , l{aj) A 2 "Q / A X 2^ = 1; . . . , ^ (3 2g) 
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Figure 3. Homology basis in the Riemann surface. 



3.4 Field equations in terms of Abelian differentials 

The characterization of the derivatives of the prepotential with respect to the partial 't Hooft 
parameters is essential to determine the prepotential from (3.12) and (3.15), as well as to 
study phase transitions of spectral curves. In Appendix A we prove that 



In this section we show that the second-order derivatives (3.28) can be expressed in terms 
of Abelian differentials of the two-sheeted Riemann surface determined by (1.5). 

We take the homology basis of cycles {A±, . . . ,A s —i,Bi, . . . ,-B s _i} as shown in fig- 
ure 3, and introduce the meromorphic differential y(z)dz, where j(z) is the extension of 
the function (1.6) to the two sheets of the Riemann surface by means of the two branches 
yi(z) = ~V2{z) = y(z), so that i/2(z-) = -y 2 (z+) = yi(z + ) = -yi(z-). Moreover, us- 
ing (3.5) and (3.8) we find that the 5-periods are 



Let us denote by d(f>i, . . . ,d^ s _i the canonical basis of normalized holomorphic differ- 
entials 




(3.27) 




(3.28) 




(3.29) 




(3.30) 



- 12 - 



We recall that these differentials are of the form 

Mz) = ^r\dz, (3.31) 
w(z) 

where Pj{z) are polynomials of degree not greater than s — 2. Likewise, we denote by d^o 
the third kind normalized meromorphic differential 

dft = 0, i = l,...,s-l, (3.32) 

whose only poles are ooi and 002, and such that 

I (—z 1 + 0{z ))dz, z — > 002. 

It can be written as 

dn (z) = ^dz, (3.34) 

where Pq(z) is a polynomial of degree s — 1. In appendix A we derive the following expression 
for the second-order derivatives of the prepotential, 

= 4vri(^ - 1)M4) - 2^o«), (3.35) 

where the corresponding Abelian integrals (f)j and f^o are defined in appendix A by integra- 
tion in the first sheet of the Riemann surface. Therefore, the field equations (1.15) admit 
the following general formulation in terms of Abelian integrals: 

s-l 

47riJ] Nj (j)j(af) + 2Nn (af) -log A 2N = 0, i = l,...,s. (3.36) 
i=i 

3.5 The case s = 1 

In the one-cut case if we denote = a, af = b, S = S\ and L = L\ then, differentiat- 
ing (3.15) with respect to 5, we find 

i = j rV W «iM + 5g, ( , 3 r, 

and 

dz 

dn (z) = (3.38) 
yj(z - a)(z - b) 



d 2 T , fb-a^ 2 



Hence (3.35) reduces to [6, 17, 18] 

Log(^— -) . (3.39) 



(3.40) 



OS 2 

Using (3.38) and (7.12) (cf. appendix A) it follows easily that 

da 4 db 4 



dS h{a){a-bY dS h{b){b-a)' 
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and then (3.39) implies 

d 3 F 8/1 1 , 



dS 3 (b-a) 2 \h(a) h{b), 

where h{z) is given by (2.5). 

Several useful results follow from these formulas. For instance, by substituting recur- 
sively (3.38) into (3.37) and into (3.15) we derive the following general expression for the 
prepotential in the one-cut case: 

t- 1 f h „r,s,,, S f b W(z)dz S 2 T (b-a\ 2 



Moreover the field equation in the case of unbroken gauge group (s = 1) 

d 2 i 

dS 2 

takes the form 



ivf5 + logA 2JV = 0, (3.43) 



log(^) =0, (3-44) 

where 5 = (b — a)/2. Thus there are N different values of 5 2 characterizing the vacuum 
states 

5 2 k = 4( k A 2 , k = l,...,N. (3.45) 

Furthermore, from (3.41) we deduce that the singular (non-analytic) solutions of the field 
equation (3.43) may only arise near points (5 C , A c ) such that the one-cut spectral curve for 
S = S c satisfies one of the conditions 

(1) The cut shrinks to a single point a = b. In this case S c = 0. 

(2) A double root of y 2 (z) collides with a cut endpoint i.e h(a) = or h(b) = 0. 

(3) It is verified that 

h(a) = -h{b). (3.46) 

4 Phase structures and critical processes 
4.1 Critical spectral curves 

In order to formulate the notion of critical spectral curves we first recall that a point of a 
real curve F(x, y) = is critical if at that point both partial derivatives dF/dx and dF/dy 
vanish. In the case of a minimal cut 7,- of a spectral curve £(7, S) we have that F = Re Gj 
and 

^-ReG 3 (z) = Re{e- i ^ s >y(z)), (4.1) 

-^-ReGj(z) = -Im(e- iarg S(z)) . (4.2) 

Therefore the cut jj has a critical point if for a certain value S c of the set of 't Hooft 
parameters a zero zq of y(z) different from meets the path jj. In this case the minimal 
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character of the cut is lost because the phase of y(z)dz at zq is undefined, and we say that 
the corresponding spectral curve £(7, S c ) is critical. As we will illustrate in the study of the 
cubic model, instances of these critical spectral curves happen in phase transition processes 
of splitting of minimal cuts. In general critical spectral curves are common limits of several 
families with different number of cuts and they arise when a zero zq of y(z) different from 
meets one of the Stokes lines emerging from a - . 

Critical spectral curves exhibit not only splitting of cuts but also birth and death of 
cuts at a distance as well as merging of two or more cuts [20]. It should be noticed that 
the merging of two minimal cuts with partial 't Hooft parameters S\ and S2 gives rise to a 
minimal cut only if arg S\ = arg S2 . 



4.2 Prepotential and its derivatives at the splitting of a cut 

Our main application of the discussion of section 3.4 concerns the behavior of the prepo- 
tential and its derivatives at a splitting of a cut. Thus let us consider a family of (s — l)-cut 
spectral curves such that the minimal cut 7 m _i splits into two minimal cuts 7 m _i and 7 m 
with 



"m-l 



a, 



(4.3) 



to give a new s-cut spectral curve. If we denote by super indices ( s ^ and ^ the respective 
magnitudes, it follows at once that 



w 



(z) = (z- a)w^{z), p^ s) {z) = p is - 1] {z). 



(4.4) 



The corresponding critical values ^ and of the partial 't Hooft parameters are 
related by 



S. 



(s-l) 



S { k s \ for 1 < k <m-2, 

S^ + Stf, for k = m-l, 
v. Sje+v fox m < k < s — 1, 



(4.5) 



where argS^_^ 



arg bm. 



In appendix B we show that the third kind 
normalized differential df^o and the normalized holomorphic differentials dcfti satisfy 



d^" 1} 



' 

.w 
k ' 



for 1 < k < s — 1, 



(4.6) 
(4.7) 



where 



k, for 1 < k < m — 2, 

m — 1 (or m) for k = m — 1, 
k + 1, for m < k < s — 1. 



(4.8) 



As a consequence of (3.10), (3.27), (3.35), (4.4), (4.6) and (4.7) we have the following 
relations for the prepotentials J 7 ^ -1 ) and J 7 ^ and their first and second order derivatives 
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at their respective critical values S^ s ^ and S^: 

jr(s-l) = jr(s)^ ( 49 ) 
aT-(s-l) <9F( s ) 

-^-^= , ferl<fe< g -l, (4.10) 

for 1 < *,i < a - 1. (4.11) 



ds\'- 1) as\'- 1) dSi s) dSi sy 

Let us consider now a family of spectral curves parametrized by a real control parameter 
T such that at a certain critical value T = T C there is a splitting of one cut, i.e for T < T c 
(T > T c ) the spectral curves have one minimal cut (two minimal cuts). Let us assume that 
at the critical value T c 

S M = sP+S?\ = s{ 2) + Si 2) , §?> = §® + S?\ (4.12) 

where dots stand for derivatives with respect to T. Then from (4.9)— (4.11) we have that at 
the critical value T c 

J-( 2 ) = JFM, (4.13) 



dS 1 



- { S[ 2) + S! 2 2) ) = — 1 (4.14) 



«(2) , V ,.,2) s ,,2) 



d 2^(2) 5JT(2) 9^(2)^ A 5^( 2 ) 



= 7^0) ^ + ^ } ) + + S 2 2) ) 2 = (4-15) 

dSj ; d(sf } ) dJ 

Thus the prepotential and its first two derivatives are continuous at S = S c . However, 
examples in random matrix theory show a jump discontinuity for the third order deriva- 
tive [10, 21, 22]. Therefore the splitting of one cut is expected to be generically a third-order 
phase transition in the space of spectral curves with minimal cuts. 

4.3 Splitting of one cut in the quantum parameter space 

For simplicity let us consider a critical point of a one-cut family of spectral curves 

(2) (2) 

corresponding to a splitting into two cuts with partial t Hooft parameters S\ and S% . 
Then from (4.11) we have that at these critical values 

d 2jr(l) d 2jr(2) 



(1K2 «q( 2 W 2 )' 



a(sr) 2 asf'ds. 



for i,j = 1,2. (4.16) 



j 



As a consequence the sectors Aiq and Adq of the quantum space of parameters touch at 
the value of A 2 given by 

A c 2 = exp( -^-4— ). (4.17) 



-16- 



Indeed, it is clear that and (S^\ S^) satisfy the field equations 

N ^S^+ l °S A c N = 0, (4.18) 

and 

2 P) 2 f( 2 ) 

i=l J 

respectively. In this way we have that a critical spectral curve corresponding to a one-cut 
splitting determines an interpolating point between the vacua spaces corresponding to the 
phases of unbroken U(N) and broken U(N\) x U(N2) gauge groups. 

5 The cubic model in the one-cut case 

In this section we apply our theoretical results of sections 2 and 3 to provide a fairly complete 
global description of the spectral curves with one minimal cut for the cubic potential, which 
we write in a form resembling the standard form of the exponent in the integral expression 
of the Airy function, 

z 3 

W(z) = — -wz. (5.1) 

O 

5.1 Endpoints, series expansions and prepotential 

We will denote = a, = b, and S = Si. Following the procedure outlined in section 2.1 
the corresponding y function to be substituted in (1.5) has the form 



y{z) =(^z+ *(a + b) J ^{z-a){z-b), 



(5.2) 



and f(z) = —ASz + &o- The resulting equations for a and b are simpler when expressed in 
terms of their semi-sum (3 = {b + a)/2 and semi-difference 5 = (b — a)/2, 

2(3 2 + 5 2 = 2w, (5.3) 

(35 2 = 25. (5.4) 

Therefore f3 satisfies the cubic equation 

/3 3 - wp + S = 0, (5.5) 

and 5 is determined by 

£ 2 = y, if P 7^0 (5.6) 

5 2 = 2w, if p = 0. (5.7) 

Note that in the latter case S must be zero. 
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Figure 4. lm/3o(S), Im/3i(5), and ImfaiS) for w — 3/2 2 / 3 . For this value of w the branch points 
are S± — ±1. 

The solutions /3(w, S) of the cubic equation (5.5) — as well as 5(w, S) and the endpoints 
a(w, S) and b(w, S) — satisfy the scaling relation 

0(w,S) = w 1/2 P(1,S/w 3/2 ), (5.8) 

so that, aside from the w 1 / 2 factor, these magnitudes are functions of the single complex 
variable S/w 3 ^ 2 . In our context it is natural to fix the value of w (i.e., to fix the potential) 
so that the solutions of the endpoint equations for the one-cut case of the cubic model are 
described by the three-sheeted genus zero Riemann surface (5.5). The corresponding three 
branches of j3 are given by 

where 

Ms) = ^^| + Vt -(!)'■ (5 - 10 > 

Here we assume that the cubic root has nonnegative real part and that the square root has 
nonnegative imaginary part. The finite branch points where two roots coalesce are 

S± = ±2(u>/3) 3/2 , (5.11) 

at which /8b(£_) = ft OS-) = and Pi(S+) = f3 2 (S+) = ^JwJ?,. 

In the three plots of figure 4 we show the imaginary parts Im /3q(S), Imf3i(S), and 
ImfaiS) respectively of the three branches of /3(S) for a fixed value w = 3/2 2 / 3 . This 
value of w has been chosen because the corresponding branch points are S± = ±1, and 
therefore easily identifiable in the figure. The upper and lower edges of the cuts along 
HeS < S- = — 1 of /3o(S) (red surface) and ^(S) (blue surface) have to be glued, as well 
as the edges of the cuts along ReS > S+ = 1 of 0i(S) (green surface) and ^(S) (blue 
surface). 
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Thus, given a value of S ^ S±, there are three values /3k(S) from which we find three 
possible pairs of endpoints, 

a k (S) = P k (S) - s/2S/MS), (5-12) 
b k (S) = P k (S) + ^2S/p k (S). (5.13) 
Using (5.9) we find the expansions of (3 k (S) as S — > 

MS) = -V^-^ + ^j- 2 +0(S% (5.14) 

A(s) = V s____ +0(n (5 . 16) 

ft(S) = -+0(S 3 ). (5.16) 
u> 

Therefore according to our discussion in section 2.3, as S — > the branches /?o and P\ 
represent families of one-cut spectral curves with a classical limit in which the cuts shrink 
to the critical point —y/w and +y/w of the cubic potential. Note also that the endpoints a 
and b can be expanded as Puiseux series. For example, 

9 l/2cl/2 c C-3/2 qc2 

«^) = ^ ± ^-^ ± w- 8 ^ + °( s5/ ' 2 )' < 5 ' 17 > 

where = a and = 6. The function ^(S) does not represents a family of spectral 
curves with a shrinking minimal cut as S — > since it reduces to the solution (5.7) with 
ai± — > ±y/2w. We will see in sections 5.3 and 6.2 that there are not spectral curves with 
a minimal cut corresponding to foiS) near S = and that the limit S — > represents a 
critical spectral curve with two minimal cuts such that Si + S*2 = 0. 

The expression for the prepotential (3.42) takes, up to an unessential constant term 
— «; 3 /12, the following simple form in terms of the function @(S) 

T ~ 12 ^ 2 L ° g J • (5 ' 18) 

This expression for the prepotential is equivalent to previous results in which the first two 
terms in the right-hand side appear in rational form. For example, equation (51) in [23] can 
be put in our polynomial form using equation (49) of [23] to eliminate the denominator. 

5.2 Superpotentials and the quantum parameter space 

The one-cut case can be used to determine the vacua structure of the unbroken gauge group 
U(N) [18]. Using (5.18) and taking into account that 

dp 1 



(5.19) 



3 v r ' ° V 2 

2 / 2 A 2 \ N 

--JV/3 3 + M-/3 3 )log(— - —j . (5.20) 



dS w- 3p 2 ' 

the superpotential can be expressed in terms of P as 

Wefr = -Inp 3 - N(wP - p 3 )Log (^r- ) +M-/3 3 )logA 2JV 
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We can obtain the vevs S since (3.45), (5.3) and (5.4) imply 

(5 (fc) )2 = 4A 2 Cfc5 p (k,±) = ±v / w _2A2 Cfe , k = l,...,N, (5.21) 

and 

s (k,±) = ± 2A 2 ( fc ^u, - 2A 2 Cfc, fc = l,...,iV. (5.22) 

Hence we find two families of iV quantum vacua \k, ±) which correspond to the two classical 
vacua i-^/u; of the cubic model. 

To calculate the low-energy superpotential W\ ow (w, A) associated with the \k, =t) vacua 
we notice that from (5.21) it follows that 

2A 2 Cfc 

1, (5.23) 



so that (5.20) leads to the following simple exact expression 



^i ( ot ±} = -f W (M:) ) 8 = T^N(w - 2A 2 Cfc ) 3 / 2 . (5.24) 



We can now calculate the vevs 



dW ( k >±) 



( Tl ty(k,±) = Jgw_ = ± N ^/ W _ 2A2( fe = NP^, (5.25) 

and check that 

Sik,±) = = 2A 2 M fc ^a. (5.26) 

In this way we have that the one-cut sector A4q^ of the quantum parameter space can 
be decomposed into N subsets 

M^=utiM^' k \ (5.27) 

(1 k) / 

where A4q ' is represented by the two-sheeted Riemann surface of the function \/w — X. 
Here 

A = 2A 2 Cfc. (5.28) 

(1 k i) 

The two sheets M\ ' ' are connected through the branch cuts emerging from A = w. Each 
point (A, ±) 6 .Mq 1 '^ 1 ^ determines a unique spectral curve given by 

S = ±\Vw - A, P = ±\/w - A, 5 2 = 2A. (5.29) 

5.3 Critical spectral curves in the complex S plane 

In this section we derive an analytic condition met by critical curves in the complex S 
plane, i.e., by the locus of values of the complex 't Hooft parameter S that yield critical 
spectral curves with minimal cuts. Let us denote by X^S") the families of one-cut spectral 
curves corresponding to the branches /3^(5) for k = 0, 1, 2. The minimal cuts of Xfc(S') for 
a fixed value of S are Stokes lines defined by 

ReG fe (z) = (5.30) 



-20- 




Figure 5. Critical curves in the sense that ReG k (-/3{S)) = for w = 3/2 2 / 3 . Each arc is marked 
with the corresponding branch j3k and drawn in the color matching the corresponding branch in 
figure 4. 



where 



G k (z) 



-i arg S 



(5.31) 



Therefore, if for a certain value of S the double zero —/3k of y(z) meets the Stokes line, i.e., 
if 

ReG fc (-/3fe) = 0, (5.32) 

the curve has a critical point. In fact, we can find an analytic condition (which, however, 
has to be solved numerically) for this critical behavior, because the integrals (5.31) can 
be evaluated in closed form. Thus, we find that the critical curve in the complex 5-plane 
corresponding to a branch ftk(S) is given by 



Re 



w 
35 



6/3| - 2w + log 




-2p k + J60l - 2w 



0. 



(5.33) 



In figure 5 we show the three branches ReGfc(— 0) = (k = 0,1,2) obtained by nu- 
merical solution of (5.33) for w = 3/2 2 / 3 . Each arc is marked and colored to match the 
corresponding branch fik in figure 4 (e.g., the red critical curve, marked /3q, lies on the first 
plot of figure 4). We have performed extensive numerical calculations of the corresponding 
Stokes graphs, from which we infer the following picture. 

Consider first the branch PziS). The critical (blue) curve in figure 5 separates the S 
plane into a bounded region and an unbounded region. In the unbounded region we always 
find a minimal cut joining a to b, and therefore a one-cut solution. The critical (blue) line 
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Figure 6. Stokes complexes on the branch (3 2 (S) for w — 3/2 2 / 3 and three values of S: S — — 0.5— i 3 
(slightly below the blue critical curve in figure 5), S = —0.5 — i 2.514 668 . . . (on the critical curve), 
and S — —0.5 — i 2 (slightly above the critical curve). 

corresponds to configurations where the double root — (32 meets the minimal cut, and in 
the bounded region of the S plane the Stokes graphs do not feature finite Stokes lines, i.e., 
there are not one-cut spectral curves with a minimal cut in this branch. This behavior is 
illustrated in figure 6, where we show the Stokes graphs for w = 3/2 2 / 3 and three values 
of S: S = —0.5 — i3, in the unbounded region of the S plane and slightly below the blue 
critical curve, S = —0.5 — i 2.514 668. . . on the critical curve, and S = —0.5 — i2 in the 
bounded region of the S plane but slightly above the critical curve. Numerical calculations 
to be presented in the next section show that these critical configurations can be continued 
to spectral curves with two minimal cuts. Thus, crossing the blue critical curve from the 
bounded to the unbounded region would correspond to a "merging of two minimal cuts" in 
the terminology of random matrix theory. 

The behavior of the branches f3o(S) and /3i(S) is the same, and qualitatively different 
from the behavior of fy^S). For concreteness we will describe the behavior of f3o(S). The 
corresponding critical curve (the red curve in figure 5) separates the S plane into a bounded 
region and an unbounded region. In the unbounded region we always find a minimal cut 
joining a to b. The critical curve corresponds to configurations where the double root — (3q 
does not meet the minimal cut joining a to b, but a second finite Stokes line joining b to 
—/3o appears. However, if we proceed to the bounded region of the S plane the Stokes 
graphs again feature a finite Stokes line joining a to b, i.e., there are one-cut solutions with 
a minimal cut. This behavior is illustrated in figure 7, where we show the Stokes graphs for 
w = 3/2 2 / 3 and three values of S: S = —0.5 — i 0.5, in the unbounded region of the S plane 
and slightly below the red critical curve, S = —0.5 — i 0.388 126 ... on the critical curve, 
and S = —0.5 — i 0.3 in the bounded region of the S plane but slightly above the critical 
curve. Numerical calculations to be presented in the next section also show that crossing 
the critical curve of the S plane in the (3q branch from the unbounded to the bounded 
region we can continue with one-cut solutions with a minimal cut (as illustrated in figure 7) 
or alternatively we can continue from the critical configuration to a spectrum curve with 
two minimal cuts. In this sense, crossing the red critical curve from the bounded to the 
unbounded region may correspond either to no phase change or to a "birth of a minimal 
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Figure 7. Stokes complexes on the branch /3q(S) for w = 3/2 2 / 3 and three values of S: S = 
—0.5 — i 0.5 (slightly below the red critical curve in figure 5), S = —0.5 — i 0.388 126 . . . (on the 
critical curve), and S — —0.5 — i 0.3 (slightly above the critical curve). 

cut at a distance" in the terminology of random matrix theory [24]. 

5.4 Spectral curves with one minimal cut in the quantum parameter space 

In the case of the cubic model each point (A, ±) 6 Aiq' k determines a spectral curve (5.29) 
which will be critical if it satisfies the condition (5.33) or, equivalently, in terms of the vari- 
able A 

w y/Aw - 6A , / y/Aw - 6A =F 2y/w - \ s 



Re 



+ log 



3A V^A V \/2A 

This equation defines a curve in M ( q A±) and since 

y/4w-6X-2y/w- X\ / ^Aw - 6A + 2^/w - A 



0. (5.34) 



2A J V V2A 



(5.35) 



f\ k 4-) (lk — ) 

the parts of the curve lying in each sheet _Mq ' ' and Ai\ ' ' are identical. Figure 8 
shows that this critical curve is composed of two ovals, one of them containing the branch 
point A = w. 

(1 k ±1 

It is not difficult to give an argument that suggests which regions of the sheets Alq ' ' 
correspond to spectral curves with one minimal cut: 



1. As A —7- then S(X) — > and /3(A) — > ^/w. Therefore, according to (5.15), the 
function /3(A) near A = behaves as (3i(S) near 5 = 0. 

2. As A — > w then S(X) —> and /3(A) — > 0. Therefore, according to (5.16), the function 
/3(A) near A = w behaves as ^(S) near 5 = 0. 

As a consequence, the analysis in sec. 5.3 permits us to conjecture that the points in the 
interior of the oval containing A = w do not supply spectral curves with one minimal cut, 
and that the oval itself describes processes of splitting of a cut. Conversely, points in the 
interior of the leftmost oval determine spectral curves with one minimal cut and points on 
the oval lead to critical spectral curves describing birth of a cut at a distance. Finally, since 
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-MR MS! MR J» fK C». CK P» 

Figure 8. Critical curve in M i q' k ' ±) for w = 3/2 2 / 3 . 

as A — > oo then S(X) — > oo so that according to the analysis in sec. 5.3 we expect that points 

outside the critical ovals determine spectral curves with one minimal cut. These arguments 

are supported by numerical calculations, two of which are illustrated in figure 9. The three 

graphs in the first row correspond to three values of A on the same vertical Re A = 0.8. 

The first value is above the left oval, the second value is critical, i.e., on this oval, and the 

third value is already in the interior of the oval. These graphs show how we proceed from 

a spectral curve with one minimal cut, through a critical configuration (in which the three 

turning points are joined by finite Stokes lines), to another spectral curve with one minimal 

cut. The three graphs in the second row correspond to three values of A on the vertical 

Re A = 2.3. The first value is above the right oval, the second value is critical, i.e., on this 

oval, and the third value is already in the interior of the oval. Now the graphs show that 

we proceed from a spectral curve with one minimal cut, through a critical configuration (in 

which the three turning points are joined by finite Stokes lines), but that as we enter the 

interior of the oval it does not exist an spectral curve with one minimal cut. 

As we discussed in section 4.3 critical spectral curves describing a splitting of a cut 

process determine interpolating points between the sectors A4q^ and A4^ of the quantum 

parameter space. Therefore for the cubic model the ovals containing A = w in the sheets 
(1 k ±) 

M\ ' ' are interpolating curves between the subsets of these sectors corresponding to 
spectral curves with minimal cuts. These results show a structure of the quantum parameter 
space drastically different from that found in [18] where, due to the absence of a concrete 
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Figure 9. Stokes complexes for two paths in Adq k,ziz ^ which enter vertically from above the left 
oval (first row) and the right oval (second row) of the critical curve of figure 8. 



choice of the cuts, the only interpolation found between the sectors Aiq^ and .Mq 2 ^ is the 
single point A = w. However, as it will be proved in section 6.2, this point describes a 
merging of two minimal cuts into a non minimal cut. 



6 Two-cut spectral curves in the cubic model 

In this section we consider the two-cut phase of spectral curves for the cubic potential (5.1). 
We first derive the series expansions around the critical points of the potential for the 
endpoints of two-cut spectral curves with classical limit. Then we discuss the particular 
solution corresponding to the slice Si + £2 = 0- Finally, we describe a numerical method 
to calculate spectral curves with two minimal cuts and their critical processes. 



6.1 Endpoints and series expansions 

In the two-cut case we will denote aj~ = a, a\ = b, = c and = d. Following again the 
procedure outlined in section 2.1 the corresponding y function to be substituted in (1.5) 
has the form 

y(z) = y/{z-a){z-b)[z-c){z-d), (6.1) 
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and f(z) = 



— 4(S*i + S 2 )z + &o- The equations for the possible endpoints turn out to be 



abc + abd + acd + bed = 4(Si + S 2 ), (6.2) 

ab + ac + ad + be + bd + cd = —2w, (6.3) 

a + b + c + d = 0, (6.4) 

I(a,b,c,d) = 2ttiSi, (6-5) 

where 

f-b 

I(a,b,c,d)= / (z + — a){z + — b){z + — c)(z + — d) dz. (6-6) 



(i 



Although (6.6) can be expressed in terms of elliptic integrals it is clear that, except in 
specially simple particular cases, the system (6.2)-(6.5) cannot be solved in closed form. 
However, we can characterize solutions that have a classical limit by means of their series 
expansions for the endpoints that tend to the critical points ±\/w of the cubic potential 
as Si and S 2 tend to zero. These series are the analog of equation (5.17) in the one-cut 
case. Since in the two-cut case we lack closed-form solutions analog to (5.9), we resort to 
direct substitution of suitable Puiseux expansion into equations (6.2)-(6.5). We do have, 
however, an scaling relation analog to (5.8): 

a(w, S u 5 2 ) = w^ 2 a(l, Si/w 3 / 2 , S 2 /w^ 2 ), (6.7) 

with similar relations for b, c and d. Therefore we can write the following Puiseux expan- 
sions: 



w 




= ^1-1 + 1^1, (6.8) 

*-^(-i+r.£&\. (^) 



1 + E^7I • («») 



where we have separated explicitly the first term to identify a and b as the endpoints of 
the cut that opens up from —y/w and c and d as the endpoints of the cut that opens up 
from yjw, and where the coefficients dj and Cj are functions of Si and S 2 - Substitution of 
these series into the first three equations (6.2)-(6.4) is straightforward. As to the fourth 
equation (6.5), which involves the integral (6.6), we can formally fix the endpoints of the 
contracting integration interval with a linear change of variable 

I(a, b, c, d) = -i(b - a) V(c -a){d-a) [ y/t(l - t)r ab (t)dt, (6.12) 

Jo 
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where 



Then we expand r a b(t) as a power series in t around t = (recall that b — a tends to zero 
while the denominators remain finite) and integrate term by term. All the integrals are of 
the form 

f 1 r-, c* AArf/c + 3/2) 

j o Vt^T)t k dt=^^ w ^l, . = 0,1,2,... (6.14) 

and the substitution of the Puiseux series for a, b, c, and d in the resulting expression is 
again straightforward. 

To illustrate the pattern of the resulting Puiseux series we show the first six coefficients 
a& and c^: 



a\ = iy25i, 

°3 = -i 2 Q Sl (2Si - 35 2 ), 

o 

a 4 = l(3S 2 l -8S l S 2 + 3Si), 

o 



a 5 = i ^r( 36S f - i225i5 2 + ms z 2 ), 

IZo 

a 6 = -^(5i - 5 2 )(165 1 2 - 595i5 2 + 165 2 2 ), 



ci = -y/2S 2 , 
c 2 = -(Si - 5 2 ), 

c 3 = ^(35i-25 2 ), 

c 4 = -^(35 2 -85i5 2 + 35 2 2 ) 

o 

c 5 = -^( 695 i 2 " 122S i^ + 365 2 ), 
c 6 = ^(5i - 5 2 )(165 2 - 595i5 2 + 165 2 2 ). 



6.15) 
6.16) 

6.17) 

6.18) 

6.19) 

6.20) 
6.21) 
6.22) 

6.23) 
6.24) 

6.25) 
6.26) 



The odd coefficients and Ck involve the square roots of the corresponding 't Hooft pa- 
rameters (and cancel in the series for the differences b — a and d — c), while the even 
coefficients are homogeneous polynomials of degree k/2 in 5i and 5 2 . Note also that if 
in the expressions for the we set S\ = and 5 2 = S we formally recover the one-cut 
expansions (5.17). 

6.2 The cubic model on the slice S\ + 5 2 = 

As an explicit application of the endpoint equations (6.2)-(6.5) we describe now the two- 
cut spectral curves of the cubic model in the Seiberg-Witten slice S\ + 5 2 = [25]. These 
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(6.29) 



configurations have been used in the literature to study the phase structure of brane/anti- 
brane systems at large -/V and non-perturbative effects in matrix models [4, 8, 26, 27]. We 
look for solutions of (6.2)-(6.5) satisfying 

c=-b, d = -a, (6.27) 

where we assume \a\ > \b\. Then these equations reduce to 

a 2 + b 2 = 2w, (6.28) 

/ yj(z\ - a 2 ) (4 -b 2 )dz = 2mSi. 

J a 

We may solve (6.28) in the form 

a = \a\e iaxsw / 2 , b = |&|e iars ™/ 2 , \b\ = y/2\w\ - |a| 2 , (6.30) 

and it is immediate to see that the straight line segments [a, b] and [c, d] are minimal cuts 
provided that 

3 

arg Si = - argw. (6.31) 

In particular (6.29) reads 

•M , 

^(x 2 - \a\ 2 )(x 2 + \a\ 2 - 2\w\)dx = 2m\Si\, (6.32) 

V2M-M 2 

which gives |<Si| as a function of \a\. 

Therefore (6.31) determines a family of spectral curves £(7, Si, S2 = —Si) with two 
minimal cuts, parameterized by the values of \a\ in the interval [y \w\, W2\w\], The limit 
\a\ — > w2\w\ of this family deserves a particular attention since the spectral curve reduces 
to the one-cut solution (5.16) determined by 

5 = 0, (3 = 0, 5 2 = 2w. (6.33) 

However due to the fact that argSi 7^ argS2 the cut resulting from the merging is not a 
minimal cut. This spectral curve corresponds to the singular point A = w of the quantum 
parameter space .Mq 1 ^ at which the solutions 5^ (A 2 ) of the field equation are not analytic, 
as it should be expected since h{z) = z + j3 and the condition (3.46) is satisfied. 

6.3 Numerical calculation of minimal cuts and critical processes 

A direct numerical solution of the system (6.2)-(6.5) without a suitable, well identified 
initial point would be extremely difficult (note that for a fixed value of w the endpoints 
are functions of the two complex variables Si and 82)- However, we can take advantage 
of our knowledge of the critical curve (5.33) and the corresponding explicit solutions for 
the one-cut endpoints given by (5.9), (5.12) and (5.13), and proceed iteratively by small 
increments in Si and S2 to calculate the solutions of (6.2)-(6.5) at any desired pair of values 
(Si, S2) using as initial approximation at each step the results of the previous one. 
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Figure 10. Stokes complexes in the two-cut region for w = 3/2 2 / 3 and two pairs (Si, S2) close to the 
critical curve: the first graph corresponds to Si = S2 = 995/200, where S = —0.5 — i 2.514 668 . . . 
is the critical value in figure 6 (note that in this case arg(Si) = arg(S f 2) and the Stokes lines do not 
cross); the second graph corresponds to Si = 995'/200 and S2 = e 17T ^ 25 Si (note that in this case 
arg(Si) 7^ arg(S l 2) and two Stokes lines cross). 

To illustrate this approach, consider the second, critical configuration in figure 6. In 
figure 10 we show the numerical continuation from this critical configuration into the two- 
cut region for two examples. In the first graph of figure 10 we have set Si = S2 = 995'/200, 
where S = —0.5 — i 2.514 668 ... is the critical value in figure 6. Note that S\ + S2 is a small 
perturbation of S (cf. equation (4.5)), i.e., we split the double zero —(3 into two simple zeros 
of y{z) 2 and the critical Stokes graph into a a graph with two finite Stokes lines of similar 
length (recall that the critical configuration is not symmetric) with arg(S'i) = arg(5 < 2). Since 
the corresponding functions Gj(z) in (2.10) have the same values of argSj, the Stokes lines 
cannot cross. Similarly, in the second graph of figure 10 we have perturbed slightly this 
solution to S± = 99<Sy200 and S2 = e 17r / 25 <Si. Now argSi 7^ argS^ and two of the Stokes 
lines corresponding to different values of j do cross. As we anticipated in the previous 
section, looking at this process backwards we have an instance of a "merging of two minimal 
cuts" in the terminology of random matrix theory. 

Likewise, in figure 11 we illustrate the numerical continuation from the critical config- 
uration in the second graph of figure 7 to the two-cut region. In the first graph of figure 11 
we have set Si = S, S2 = S/10, where S = —0.5 — i 0.388 126 ... is the critical value in 
figure 7. Again, this is a small perturbation of the critical configuration which we use to 
generate the initial approximation. The double zero — (5 splits into two simple zeros and 
the critical Stokes graph into a graph with a long finite Stokes line corresponding to Si 
and a short finite Stokes line corresponding to S2. Since in this first graph we have taken 
arg(Si) = arg(S"2), the corresponding functions Gj(z) in (2.10) have the same values of 
argSj, and the Stokes lines cannot cross. In the second graph we have increased the argu- 
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Figure 11. Stokes complexes in the two-cut region for w = 3/2 2//3 and two pairs (Si,^) close to 
the critical curve: the first graph corresponds to Si — 5, S2 — 5/10, where 5 = —0.5 — i 0.388 126 . . . 
is the critical value in figure 7; (note that in this case arg(5i) = arg^) and the Stokes lines do 
not cross); the second graph corresponds to 5j = 5 and S2 = e 17r / 30 5/10 (note that in this case 
arg(Si) 7^ arg(52) and two Stokes lines cross). 

merit of S2 by -/r/30, and two Stokes lines corresponding to different values of j cross. As 
we anticipated in the previous section this is an instance of a "birth of a minimal cut at a 
distance" in the terminology of random matrix theory. 

Therefore we have an efficient numerical algorithm that taking critical one-cut solutions 
as data to generate an initial approximation to solve equations (6.2)-(6.5), allows us to 
proceed stepwise tracking the solution along any specified path of the complex variables 
Si and 52, generate the two-cut endpoints o(5i,52), &(5i, S2)) c(5i,52), d(Si,S2) and 
calculate the corresponding Stokes graph. 

7 Summary 

In this paper we have presented a self-contained approach to study the phase transitions 
that arise in the study of large iV dualities between jV = 1 SUSY U(N) gauge theories and 
string models on local Calabi-Yau manifolds. Conceptually, it is based on the use of spectral 
curves with branch cuts that are projections of minimal supersymmetric cycles. From the 
practical point of view, the key elements of our approach are a system of equations for the 
branch points and a characterization of the branch cuts as Stokes lines of a suitable set of 
polynomials. The system of equations for the branch points is derived from the form of the 
equation of the spectral curve with a fixed number of cuts after imposing period relations on 
the cuts. In turn, these branch cuts are a natural generalization of the cuts that define the 
eigenvalue support in holomorphic matrix models (note that in this latter case the partial 
't Hooft parameters are essentially the eigenvalue fractions and therefore real and positive) . 
However, we do not rely on any underlying matrix model but characterize the (in general 
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complex) density as an extremal of a prepotential naturally associated to the spectral curve. 
By writing the derivatives of the prepotential with respect to the 't Hooft parameters in 
terms of Abelian integrals we study the splitting of a cut in the space of spectral cuts and 
show it is typically a third order transition. We also show that a critical spectral curve 
corresponding to the splitting of a cut represents an interpolating points between different 
sectors in the quantum parameter space. We have applied these theoretical results and 
numerical calculations to study the cubic model, finding the analytic condition satisfied by 
critical one-cut spectral curves, and characterizing the transition curves between the one- 
cut and two-cut phases both in the space of spectral curves and in the quantum parameter 
space. Although our analytic results in the one-cut phase of the cubic model rely on the 
explicit solution of a cubic equation, in principle nothing prevents the use of the same 
numerical methods used in the two-cut phase to study spectral curves for non polynomial 
potentials. 
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Appendix A 

In this appendix we collect the proofs of the main identities satisfied by the prepotential 
associated to spectral curves. To prove (3.27) and (3.28) we differentiate (3.10) with respect 
Si taking into account that p(af) = ■ ■ ■ = p(af) = 0. Thus we have 




(7.1) 



where we have used (3.8) and the fact that because of (1.10) 




(7.2) 



Moreover, from (3.8) we have that 




7 



(7.3) 



and (3.27) and (3.28) follow immediately. 
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We will next prove that the asymptotics of Lj in the classical limit is 

Lj ~ W{aj) ~Y.Sk log A% + Sj (l + log (^^) ) , (7.4) 



where Ajk = aj — a^. To this aim we set 2j = a^" in (7.3) and get 



3 ' 2vri 
The classical limit of the first term is 



-j = W(a+) - -L I Log(z - a\) 2 y{z + )dz. (7.5) 



W"(a 



W(a+) = W(Pj + 6j) ~ W( aj + ~ W( 0i ) + — (7.6) 
The second term decomposes as a sum 

f Log(* - a+) 2 y(z + )dz = £ / Log(z - a+) 2 y(z + )dz. (7.7) 



k J ^ 

Taking into account the periods (1.7), for k ^ j we get 



f Log(z - a+) 2 y(z + )dz ~ Log(A 2 fc ) f y{z + )dz = 2vriS fc Log(A 2 fc ), k ± j. (7.8) 



Moreover, using (2.16) we obtain 



f hog{z - a+) 2 y{z + )dz ~ iW"( % ) f ' Log(z - a+) 2 J(z - a 7)(a+ - z)dz. (7.9) 
y 7j J & . 

The integral in the right-hand side can be exactly calculated using the change of variables 
z = aj + (at —aj)t: 

{ Log(z - a+) 2 J{z-aj)(a+-z)dz 

J a 

j 

= 4o 2 Log(4o 2 ) I y/t(l - t)dt + 4o 2 [ Log(t - 1)V*(! ~ *) d * 
Jo Jo 

5 2 5 2 
= tt^ Log(4o 2 ) + tt^-(1 - Log 16) 

5 2 5 2 

= 7r^(l + Log£) (7.10) 

Hence using (2.19) the asymptotics (7.4) follows. 

Let us prove now the expression (3.35) of the second derivatives of the prepotential in 
terms of Abelian integrals. Differentiating (7.3) we find that 

~ dU - W"^) 2 ^, Zi^i, (7.11) 



dSidSj dSj J oy " dS 3 
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where the integrals are independent of the choice of Z{ in j{. From (3.5) and (1.7) it follows 
that 

^-y(z)dz = -4vri(l - 8 ja )&<j>3 - 2dft . (7.12) 

Hence 

d M z l = _ 2 (l _ 5 JS )d^(z + ) - -<m (z + ). (7.13) 

Obj 7T1 

Substituting (7.13) into (7.11) we get 

d 2 T 



dSidSj 



2(1 -S js ) [ Log(z-z l ) 2 dcP j (z + ) 

+ - / Log(z - Zi ) 2 dn (z + ), Zi E 7i . (7.14) 
Consider now the first integral in the right-hand side of (7.14) and denote 

Fj(u)= Log(z -u) 2 dcf) j (z + ), UE7. (7.15) 
It follows from (3.31) that 

d</>j(z+) = -d(j)j{zJ), z G 7, (7.16) 

which implies that 

Fj-(u) = / (log(z+ - u) + log(z_ - u))d^j(z + ) 

^ log(*-u)<%(*), (7.17) 



where ^4 = + • • • + is the sum of the contours in figure 1 . Therefore 

= f ^M = , u 6 7, (7.18) 

because the integrand is analytic outside A and has residue zero at oo. Therefore the 
functions Fj(u) are constant on any connected piece of 7. Moreover, equation (7.17) shows 
that the functions Fj (u) are analytic in C \ T and 

Fj(!l , = / d Mi) = -2 I iWM, uij. (7.19) 
J J A z-u w{u) 

Since 

/ = 0, (7.20) 

it also follows from (7.17) that Fj{u) = 0(l/u) as«4 00. Therefore, using the standard 
definition of the Abelian integrals 



ru 

4>j{u) = I d(f)j(i 
J 001 



(7.21) 



-33- 



and equation (7.19), we may write 

Fj(u) = -2m<f>j(u), ueC\r, (7.22) 
and by continuity we conclude that 

Fj(zi) = -2vri^(a+), z { £ 7i . (7.23) 
The analysis of the second term in equation (7.14) is similar. We denote 

F(u) = I Log{z-u) 2 dn (z + ), u£ 7 . (7.24) 
From (3.34) we deduce that 

dtt (z + ) = -dQo(z-), 267, (7.25) 

and the same arguments used for Fj(u) show that F{u) is constant on any connected piece 
of 7, and that the derivative F'(u) outside 7 is 

F'{u) = -2vr i u d 7 . (7.26) 
w(u) 

Taking into account that 

I dn (z) = -7ri, (7.27) 
if we define the the Abelian integral Qq (it) as 

no(«)= lim ( f dfio(z) + logu ) , (7.28) 

and use equation (7.26), we find 

F(u) = -2vrif} (u), !i£C\r. (7.29) 
Therefore, by continuity, we finally get the result 

F{zi) = -2mn (at), Vzj 6 7f . (7.30) 
Equation (3.35) is now a direct consequence of (7.14), (7.23) and (7.30). 

Appendix B 

The identities (4.6)-(4.8) which describe the behavior of the differentials df^o and d0i, . . . , d0 s _i 
under coalescence of the endpoints of two cuts can be proved using a method due to 
Tian [28]. For conciseness we provide the proof for dfio only This differential can be 
written as 

dn (z) = I ^4dz, (7.31) 
w(z) 
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where Pq{z) is a polynomial of degree s — 1 which depends on the endpoints = 
(of, . . . , af) and is uniquely characterized by (3.32) and (3.33). Let us first prove that (4.3) 
implies 

0. (7.32) 



p (s W*); 



From (3.32) we have 



so that 



^ a(S) ) d , = 0, (7.33) 



a ™-i P (s) (z;a( s )) 



u;( s )(z;a( s )) 



P W(z;aW; 



d2 



(«)+ (<o- 

a„_i=a™ =a 



* m — 1 "m 



-dz, 



(z-a)^^- 1 )^^^- 1 )) 
which implies (7.32). Thus, for a coalescence of the type (4.3) the function 

P W (z;aW) 



z — a 



(7.34) 



determines a polynomial Pq S (z; a^ s 1 - ) ) of degree s — 2 and dfig reduces to 



d^(,) = 44^d*. (7.35) 
This expression determines the differential dfig* ^(z) and therefore (4.6) follows. 
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